A new class of quantum bound states: diprotons in extreme magnetic fields 
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This paper considers the possibility that two charged particles with an attractive short-ranged 
potential between them which is not strong enough to form bound states in free space, may bind in 
uniform magnetic fields. It is shown that in the formal limit where Coulomb repulsion is negligible 
(q — > and Bo — > oo with qBo fixed, where q is the charge and Bo the field strength) there always 
exists a bound state for a system of two identical charged particles in a constant magnetic field, 
provided that there is a short-range uniformly attractive potential between them. Moreover, it is 
shown that in this limit any potential with an attractive s-wave scattering length will posses bound 
states provided that the range of the potential is much smaller than the characteristic magnetic 
length, r = {^)~ 1/2 . For this case, the binding is computed numerically. We estimate the size of 
the magnetic field needed to approximately reach a regime where the formal limit considered becomes 
a good approximation to the dynamics. These numerical estimates indicate that two protons in an 
extremely strong magnetic field such as might be found in a magnetar will bind to form a diproton. 
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I. INTRODUCTION 

In nonrelativistic quantum mechanics, uniformly at- 
tractive potentials in one or two dimensions always allow 
the existence of bound states. While this is not gener- 
ally true in three dimensions, the presence of a magnetic 
field can alter the situation considerably: attractive po- 
tentials between two charged particles with insufficient 
strength to bind in free space can bind in the presence 
of a strong external magnetic field. This novel class of 
quantum bound states is the focus of the present paper. 

Consider the nonrelativistic quantum dynamics of two 
charged particles in a uniform magnetic field. It is useful 
to specialize to the case where the charge and mass of the 
two particles are identical; in this situation the center-of- 
mass motion decouples from the relative motion. We 
begin the analysis working in a formal limit: 



q — > and Bq — > oo with qBo fixed. 



(1.1) 



This limit has the property that the interaction between 
each particle and the external magnetic remains of fixed 
strength while the Coulomb repulsion between the parti- 
cles is suppressed. 

As a first example we consider a finite-ranged central- 
force potential between the two particles which is attrac- 
tive at all distances but is not strong enough to bind 
the particles in the absence of the magnetic field. In 
such circumstances there is a general theorem that the 
two particles will necessarily have a bound state even 
if the potential is not strong enough to bind in free 
space. However, the emergence of bound states in mag- 
netic fields from potentials not strong enough to bind in 
free space is more general than this. One important case 



where such binding is guaranteed to happen in the limit 
q — ► and Bq — > oo with qBo fixed is when there is a fi- 
nite ranged interaction with both attractive and repulsive 
regions but which has net attraction at low energies in 
the sense that the scattering length, a, is attractive and 
the characteristic magnetic length scale (^j 2 -) -1 ^ 2 (where 
we work in natural units natural units of % = c = 1) is 
much larger than the range of the interaction R. 

There are simple intuitive reasons why such binding 
might be expected. Consider first classical scattering in 
the situation described. Qualitatively, if the Coulomb 
repulsion is neglected, two non-interacting particles with 
the same charge and mass moving in a magnetic field will 
have the same angular frequency. Thus after they scatter 
they will bend in the same manner and thereby remain 
close to one another for a longer time than an analogous 
scattering in free space. This in turn allows the attrac- 
tive but finite-ranged interaction to remain significant for 
a longer period leading to increased net attraction. The 
quantum analog of this is presumably also increased at- 
traction compared to the free space case and it is not 
surprising that this can bring about binding. 

It can also be understood qualitatively from a quantum 
mechanical perspective. Recall that a charged particle in 
a constant magnetic field of strength Bq moves at the 
cyclotron frequency lo c — ^—^ with an energy spectrum 
described by the well-known Landau levels. This can be 
easily seen via an appropriate choice of gauge (the sym- 
metric gauge) in which case the Hamiltonian of a single 
charged particle in a magnetic field can be written as the 
sum of a two-dimensional harmonic oscillator Hamilto- 
nian and an angular momentum term in the direction 
of the magnetic field. In this gauge, the magnetic field 
(which without loss of generality we take to be in the z 
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direction) acts to restrict the particles to confined har- 
monic motion in the x — y plane, leaving free particle 
motion only along the direction of the magnetic field. 
Under these circumstances, it is highly plausible that a 
situation similar to the one dimensional case might arise 
in which any purely attractive potential between the two 
particles will generically bind them. 

Given the generality of the result, one might think that 
such bound states should be ubiquitous. However, there 
is a fundamental limitation: the analysis was carried out 
in the regime q — > and -Bo — *■ oo with qBo fixed. In 
practice, however, one cannot find particles with arbitrar- 
ily small charges since the charge is quantized in units of 
e. Thus to approach the regime of relevance one needs to 
find situations in which the Coulomb interaction in the 
lowest lying state is very small compared to the charac- 
teristic energies in the problem arising from the particles 
interacting with the magnetic field, (i.e. the cyclotron 
frequency — ui c = ^-p) and is also much smaller than the 
characteristic energy scale of the short-ranged attraction 
between the particles. In practice, these conditions are 
hard to obtain. One requires very large magnetic fields. 

The requirement of extreme magnetic fields implies 
that the effect in its simplest form will not be directly 
relevant for systems of atomic or molecular ions. Note 
that, although one can describe the interactions between 
such ions quantum mechanically via an effective poten- 
tial, the potential depends on the internal structure of the 
ions. The internal structure will be greatly altered due 
polarization effects in very strong fields. Thus one cannot 
describe the interaction of the two ions in the presence 
of the magnetic field as essentially the free space inter- 
action. Of course, there may well be situations in which 
two atomic or molecular ions bind in the presence of a 
magnetic field which do not bind in free space. However, 
such an effect would be a many-body effect in which the 
effective two-body potential between the atoms is altered 
due to the changed internal structure. This is not the ef- 
fect discussed here. As we will show, the magnetic fields 
needed to see the phenomenon discussed here are many 
orders of magnitude larger than the characteristic field 
strength needed to alter the interparticle potential by al- 
tering the internal structure 

Thus, the most natural place to look for the effect 
is where strong interactions dominate — i.e. in nuclear 
physics. The most favorable case we know is for bound 
states of the two protons. As we will discuss below, 
such diprotons are expected to form in magnetic fields of 
10 12 T. Such fields are found in nature but static fields 
of this scale are only seen in astrophysical contexts such 
as magnctars. It remains an open question as to whether 
there are any astrophysical consequences of these states 
which have direct implication for observation. Whether 
this turns out to be the case or not, these states are inter- 
esting quantum mechanical curiosities in their own right. 

Before proceeding it is useful to define what we mean 
by a bound state in this context. In free space this is 
unambiguous: the relative and center of mass motion 



decouple, and bound states have a discrete energy spec- 
trum, negative energy and a normalizable wave function 
for the relative motion; the probability of finding the par- 
ticles at some fixed distance drops off exponentially with 
distance. In contrast, continuum states in free space 
have a continuum spectrum, positive energy and non- 
normalizable states. In a magnetic field things arc dra- 
matically different. In the first place, relative and ccnter- 
of-mass motion do not generally decouple. Secondly, the 
distinction between discrete and continuum spectra be- 
comes ambiguous: even non-interacting single particles 
have a spectrum with discrete Landau levels (with infi- 
nite degeneracy) along with a continuum index. Given 
this situation, we begin as a technical matter with a def- 
inition of a "bound state" in terms of the behavior of the 
drop off of the probability amplitude with distance: a 
bound state is an energy eigenstate in which the probabil- 
ity amplitude for finding two particles a distance r apart 
falls of at long distances faster than a power law in all 
directions. This definition is, of course, fully consistent 
with zero-field case and is a useful generalization. 

In this paper, we restrict our attention to situations 
in which the charge and mass of the two particles are 
identical. As will be discussed below, provided that the 
dynamics are non-relativistic, the center-of-mass and the 
relative motion cleanly separate in this case. The center- 
of-motion behaves as a single particle of mass M = 2m 
and charge Q = 2q; it has a spectrum of (infinitely degen- 
erate) discrete Landau levels plus a continuum associated 
with the motion in the z direction. Bound states in this 
circumstance have two properties analogous to the free 
space case: i) the wave functions for the relative motions 
are normalizable and ii) the energy of the state is nega- 
tive by a non-zero amount relative to what it would be 
if the particles were noninteracting or equivalently were 
very far apart. 

The goal of this paper is to explore the phenomenon 
of binding due to intense magnetic fields rather gener- 
ally. We will show that in certain circumstances the phe- 
nomenon must occur. However, for any particular inter- 
action whether or not binding occurs and, if it does, the 
value of the binding energy requires a detailed knowl- 
edge of the interaction and a numerical solution of the 
Schrodinger equation. It is useful to try to understand 
how this effect occurs rather generally without relying on 
too many details of the interaction (which in many cases 
is not known with precision in any event). Accordingly, 
we focus here on some limiting cases through which one 
can gain insight into the phenomenon without specifying 
the details of the particle-particle interaction. The key 
limitations in this analysis are a restriction to situations 
in which the Coulomb interaction may be neglected and 
where the magnetic fields are small enough so that the 
characteristic magnetic length scale is much longer than 
the range of the potential. While these limitations are 
severe, there is at least one interesting physical system — 
a bound state of protons in magnctar strength fields — for 
which the approximations appear to be quite good. 
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The plan of the paper is as follows: We begin by set- 
ting up the general problem and then prove that, in a 
constant magnetic field and in the limit q — > and Bq — > 
oo with qB fixed, two charged particles with the same 
charge and mass and with a uniformly attractive poten- 
tial between them always have a bound state. We also 
prove that in this regime any potential which is attrac- 
tive in the sense of having an attractive scattering length 
must have a bound state in a sufficiently weak magnetic 
field and derive an analytic expression for the binding en- 
ergy for weak magnetic fields. In connection to this weak 
field limit we note that it is very difficult to simultane- 
ously have sufficiently weak fields to be in this regime 
and to have the Coulomb repulsion be small enough to 
neglect. 

Next we consider systems for which the scattering 
length is "unnatural" , ie. much larger then the char- 
acteristic range of the interaction. In this regime there 
is an effective field theory [j], Q describing the system 
and one can compute binding without a detailed micro- 
scopic knowledge of the interaction. This regime is of 
phenomenological interest in low energy nuclear physics 
since the nucleon-nucleon scattering length is much larger 
than the characteristic range of the nuclear force. We nu- 
merically calculate the binding energy of the bound state 
in this effective theory as a function of a dimcnsionlcss 
parameter proportional to aB\J 2 where a is the scatter- 
ing length and Bq is the strength of the magnetic field. 
We discuss corrections to this result both due to Coulomb 
repulsion and to short-distance effects. For the case of 
two protons we use dimensional arguments to argue that 
a regime exists in which both of these corrections are 
simultaneously moderately small. 

Following this we discuss the case where the two parti- 
cles have spin. There is a general argument that energy 
of interaction of the spins with the magnetic field are of 
the same order as the binding effects considered here. It 
is easy to sec that in the regime of validity of the effective 
theory the spin degrees of freedom decouple from the spa- 
tial ones. This means that the bound states considered 
here will continue to exist. However, in the case of iden- 
tical fermions it may well happen that even when there 
is a bound state of the sort discussed in this paper, the 
ground state of the system need not be bound. This can 
come about because the Pauli principle combined with 
the dynamics of the bound state can require the spins to 
be in distinct m-states. The energy gained in binding can 
be less than the energy lost by aligning the spins in the 
field. This situation is radically different from what oc- 
curs in free space where the ground state must be bound 
if any bound states exist. As we will see later, the situa- 
tion where the lowest bound state is not the ground state 
is precisely the circumstance for the diproton. 

We end the paper with a brief concluding section in- 
cluding a few words about the possible astrophysical sig- 
nificance. 



II. BOUND STATES IN MAGNETIC FIELDS 
A. Generalities 

In this section we set up the general problem and prove 
the existence of a bound state when there is a poten- 
tial which is attractive everywhere between two identical 
charged particles in a constant magnetic field in the limit 
q — > and Bq — > oo with qB^ fixed. 

Consider the case of two charged nonrelativistic par- 
ticles (with the same mass m and the same charge q) 
interacting via a central potential. For simplicity we will 
neglect spin effects here and return to the question of spin 
later in the paper; accordingly we will assume a central 
potential initially. For simplicity we will assume that 
the two particles have identical charges and masses. The 
Hamiltonian of the system is then 

- (pi-qA) 2 (p 2 -qA) 2 

H = 1 h V (|ri - r 2 |) . (2.1) 



2 m 



2m 



Without loss of generality we can take the magnetic 
field to be in the z direction. We find it convenient to 
work with the symmetric gauge 



yx) . 



(2.2) 



For the case of interest here where the mass and charge 
of the two particles are equal, the center-of-mass and 
relative motion separate. With fj, = ^ as the reduced 
mass and M = 2m as the total mass, the relative and 
center-of-mass Hamiltonians in configuration space are: 



Hrel = ~^V 2 + \^ 2 (* 2 + y 2 ) - uL z + V{\r\) (2.3) 
HcoM = -^"V 2 + X -Muj 2 (x 2 + y 2 ) - ljL z , (2.4) 



where u> 



qBp _ qBg 

2 l ~ c M 4/i ■ 

Iii the absence of a short-range potential, the center 
of mass and the relative equations closely resemble the 
Hamiltonian of a single particle with a Landau level spec- 
trum. In the absence of interactions between the particles 
the relative motion does as well. Therefore, without the 
potential, the total energy for the ground state of the 
relative motion is u> — A bound state corresponds 

to an energy below this value. 

It is easy to convince one's self that to minimize the 
energy the term proportional to L z in Eq. (|2.3|) should 
be zero. It is always possible to do this since L z com- 
mutes with the Hamiltonian. With the choice of gauge 
of Eq. (|2.2p the Hamiltonian can be written: 



2/i 



f 



2 2 
flUJ p 



V(f) 



(2.5) 



where V(r) is the potential, and V 2 is the three dimen- 
sional Laplacian. A two-dimensional harmonic oscillator 
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term, H2D-HO-, may be isolated from Eq. (|2.5p and the 
Hamiltonian may be rewritten as follows: 

H rel =H 2D - HO -i^-f^-V(r). (2.6) 



B. A theorem 

A variational argument may then be used to show that 
a uniformly attractive potential always binds with the 
Hamiltonian in Eq. (|2.6p . Before doing this note that 
the condition in Eq. is implicit in taking a uniformly 
attractive interaction; there will always be a long ranged 
repulsive Coulomb potential away from this limit. To 
proceed we take a simple a wave function ansiitz of 

W) = \V)2D-H0\<t>) Z , (2-7) 

where \0)2D-ho is the ground state of the two- 
dimensional harmonic oscillator in the x — y plane and 
\4>)z is an arbitrary state in the Hilbert space for motion 
in the z direction with <fi(z) = (z\4>) z ; an outer product 
between is implied. The expectation value of the energy 
in this state is: 

(i>\H rel \i>) 

= u> + J dz<f>*(z)(-^-^-V cS (zmz), (2.8) 
with V e s(z) = J dxdytpo(x,y)(V (f))ipo(x,y) 

where u is the ground state energy of the two-dimensional 
harmonic oscillator and ipo(x,y) its normalized wave 
function. Note that V'oV'o is non-negative; thus, if V(r) is 
non-negative for all r then V e s is non-negative for all z. 
That is, a purely attractive V implies a purely attractive 
V eS . 

Note that the integral in z is the expectation value for 
the energy of a one-dimensional system with potential 
Vcff. The effective one-dimensional system must have at 
least one normalized bound state with negative energy 
since there is a theorem that all purely attractive one 
dimensional systems do Thus, we can choose <f>(z) 
to be the normalized ground state of this system. In 
that case we have {il>\H re i\tp) = to — B^f where B(?, the 
binding energy of the effective one-dimensional problem, 
is positive. Thus we have found a variational state with 
energy less than that of two non-interacting particles in 
the magnetic field; the true state must be lower still and 
thus a bound state must exist. 



III. THE WEAK FIELD LIMIT 

The theorem in the preceding section is of interest as it 
illustrates the general phenomenon of particles which do 
not bind in free space forming bound states in magnetic 



fields. It is limited in that it does not tell us whether a 
potential which has both attractive and repulsive regions 
binds in a magnetic field. As a general rule one cannot 
determine whether this occurs for a particular potential 
and a given field strength without solving the Schrodingcr 
equation explicitly. 

However, there is one general situation where the exis- 
tence of bound states can be determined by general argu- 
ments: the case in which the natural length scale associ- 
ated with Landau orbits in the magnetic field, rn = i— , 

is much larger than the range of the potential between 
the particles (which we will denote as R). In this weak 
field situation, it is easy to show that bound states exist 
providing the scattering length is attractive and Coulomb 
repulsion is negligible. Moreover, if the scattering length, 
a is "natural" (i.e. of order R) and thus much smaller 
than ro, one can directly compute the binding energy. It 
is given by: 

Eb = 2(j£j u. (3.1) 

We note that the weak field condition ro R can 
be at odds with the regime in which the formal limit 
q — > and Bo — ► oo with qBo fixed (in which Coulomb 
repulsion may be neglected) can be regarded as approxi- 
mately satisfied. The neglect of Coulomb requires Bo to 
be sufficiently large; the condition ?*o R requires Bq 
to be small. For the problem of the diproton there does 
exist a window in B for which the system is in both 
regimes. 

The weak field regime is of interest in that the system 
is insensitive to the detail of the potential — at leading 
order, all that is relevant is the scattering length. This is 
precisely the situation for which an effective field theory 
[l], 0] is useful. We begin by considering the "natural 
scattering length" regime of a ~ R. We consider the 
unnatural case a 3> R in a subsequent section. 

Our goal is to demonstrate the relationship between 
the binding energy and the scattering length of Eq. (|3.ip . 
Wc begin with a general expansion of the bound state in 
a convenient basis: 

\<P) = J2\n) 2D - HO \<f>n)z (3.2) 

n 

where as before an outer product is implied; \4> n )z rep- 
resents the particular z space wave function coupled to 
a given state of the two dimensional oscillator. The key 
to our demonstration here is that in the limit a, R -C ro 
only the n = term contributes. We represent the total 
energy of the system as u>, the energy of the two dimen- 
sional harmonic oscillator, minus the quantity of inter- 
est: H re i\ip) = (— Eb + u>) \tp). Using the expansion of 
Eq. (j3~2")) and multiplying Eq. (|2~B)) on the left by the 



■5 



two dimensional harmonic oscillator state (m(p)\ yields: 
fid 2 \ 

= (-E b - nuj) (j)n{z)5 mn , (3.3) 



where V mn (z) = J dp2np tp* H (p)V(p, z)ip n (p). 

Note that the solution to the Schrodingcr equa- 
tion along the z-axis decays exponentially according to 
e -V'2t*E b z ag ^ ne b ma ; m g becomes small. Clearly, the 
n = state is the only one which contributes in the weak 
binding limit because of the (—E b — nuj) term; for nonzero 
n, the wavefunction no longer decreases exponentially in 
the i direction as the binding energy becomes very weak. 
The only term remaining that connects orthogonal m and 
n states is V mn . However, the range of the attractive 
potential V is small compared to that of the harmonic 
oscillator. Therefore, all harmonic oscillator states ap- 
pear essentially flat in the range of the potential, and all 
contribute approximately equally. 

Note that there are now three length scales in the prob- 
lem: the range of the attractive potential, R; the length 
scale of the harmonic oscillator, ?'n — —)=\ and the scale 

of exponential decay of the wave function in the z direc- 
tion, = (2pE b ) 
needed for Eq. (|3.6f 



2 / 2 . The hierarchy of length scales 
to be valid is 



R < r < k~ 



(3.4) 



One must test ex post facto to see whether this condition 
is satisfied. Note that if Eq. (|3.1|) is correct, then the 
condition in Eq. (|3 .4[) becomes 



R < r < r - 



1 



a^/puj 



(3.5) 



Eq. (J373J) now takes the form of a one-dimensional 
Schrodingcr equation. Moreover, at small binding the 
spatial extent of the one dimensional wave function 
means that the system can be well represented as a sys- 
tem with a Dirac delta potential in z with strength a and 
energy -E b : 



2/i dz 



a5(z) ) \<j) (z)) 



with Eh 



pa 
~2 



Eb\M*))> 

2 



(3.6) 



Clearly, to proceed we must determine a in terms of the 
scattering length. 

We know a priori that a does not depend on the details 
of the short distance potential but only on the scattering 
length. This suggests a useful strategy for computing a 
in general, namely to consider any class of potential for 
which we know how to compute both a and the scattering 
length and then relate the two. There is a straightforward 



situation for which we can do this: namely the case of 
arbitrarily weak potentials. For the bound state problem 
in a magnetic field this regime gives 

a = J dz2n pdpV(p,z) ^(0)ip o (0) (3.7) 

For the scattering, the weak potential limit justifies the 
the Born approximation. The Born approximation for 
the scattering length is 



a = 4irp / d xV 



(3.8) 



Combining Eqs. $3Jj§, (|3J]) and j3U) yields Eq. d3TTj) . 

A couple of quick comments about this result are in 
order. The first is that Eqs. ([3~7]) and ([3~8]) only hold 
for weak potentials. Never-the-less the result in Eq. (|3.ip 
is fully general for R,a<^ ?'o . The reason is clear — 
we have an effective theory in which the only relevant 
information about the potential is through the scattering 
length. How the scattering length arises microscopically 
is irrelevant. The second is that this result generalizes the 
theorem considered above. Any potential with attractive 
scattering length will form a bound state in a magnetic 
field provided that the system is in a regime in which 
the Coulomb repulsion is negligible and the range of the 
potential and its scattering length are much smaller than 

It is important to note however, that this weak cou- 
pling limit is of more theoretical then practical signif- 
icance. Recall that the analysis depended on two as- 
sumptions. The first was that the Coulomb repulsion 
was negligible and the second that the magnetic length 
scale is much longer then all scales associated with the 
potential including the scattering length. Clearly neither 
assumption is exact and hence there must be corrections 
to these. One can estimate the size of these corrections. 

The characteristic distance between the particles is 
(pEh^ 1 / 2 since the system is weakly bound and the wave 
function is largely beyond the range of the potential. 
From this one estimates the characteristic scale of the 
perturbative Coulomb correction to be 



/\ ^Coulomb 

E~ b 



O 



q 2 fi r'l 



(3.9) 



where q is the charge and in the second form we are 
taking the scattering length to be of the same order as 
the range of the potential as one expects for a "natural" 
(i.e. not fine tuned) problem. It is nontrivial to find a 
regime for which this is small enough to be considered 
to be a small correction. In the first place note that the 
second expression for the ratio of the correction to the 
dominant part is proportional to q 2 pR. For two light 
nuclei and strong interactions q 2 pR ^ §x- I n compar- 
ison, in atomic systems apR is many orders of magni- 
tude larger. As noted earlier this implies that binding 
of the sort considered in this will not occur in atomic or 



6 



molecular ions. In fact, it is typically marginal at best 
for nuclear systems in the weak coupling regime. Note 
that Eq. (|3.9p contains a factor of r^/R 2 . Recall that 
the weak field regime has Tq/R 2 3> 1. Thus, there is at 
best a very small window, and more likely no window at 
all, for which the system has a natural scattering length 
and simultaneously is at weak field and has perturba- 
tive Coulomb corrections. Fortunately, there are systems 
with unnatural scattering lengths where binding due to 
the magnetic field can occur. 



IV. UNNATURAL SCATTERING LENGTHS 

In the previous section we considered binding in the 
regime where both the range of the potential and the 
scattering length are much smaller than the characteris- 
tic magnetic length scale. In this section, a more general 
case is considered: one in which as before r$ 3> R so 
that the detailed structure of the short-ranged potential 
remains unimportant but where we no longer require the 
scattering length to be much shorter than rQ. Of course, 
one generally expects that in the absence of fine tun- 
ing the scattering length should be of the same scale as 
the range of the potential (or less in the case of weak 
potentials). Under such circumstances rg is also much 
larger than a. However, there are systems for which 
the scattering length is "unnaturally" large, i.e. much 
larger than the range of the potential. The case of un- 
natural scattering lengths is critical in nuclear physics 
since the nucleon-nucleon scattering lengths are unnatu- 
ral: a ~ 20 fm while R ~ 2fm Indeed, much of the 
development of effective field theory methods in nuclear 
physics in the past decade has concerned the treatment of 
systems which are unnatural [f| 0, 0, • It is important 
as a practical matter to consider the regime of unnatural 
scattering lengths. As noted above the regime considered 
previously with ro a turns out to be of very limited 
utility. 

To proceed it is useful to reexpress the problem by 
adding and subtracting a one dimensional harmonic 
potential in the z direction from the Hamiltonian in 
Eq. (|2.3p . The Hamiltonian can be rewritten as a three 
dimensional harmonic oscillator Hamiltonian minus a one 
dimensional harmonic oscillator potential term in the z 
direction: 



Hrrl — V 

2/i 



1 2 2 



\^ 2 z 2 ~V{\r\) . (4.1) 



This formulation is convenient for numerical studies in 
that one can easily implement the calculation in a trun- 
cated three dimensional harmonic oscillator basis. 

As a first step towards considering the general case 
with ro > i? but with the ratio of ro/a arbitrary, we 
consider the extreme case in which the short range poten- 
tial in the absence of the magnetic field has a scattering 
length sufficiently long to be treated as though it were 
infinite scattering (i.e. as though the system had a zero 



energy bound state corresponding to the boundary be- 
tween having bound and unbound states in the absence 
of a magnetic field). This regime is of interest theoreti- 
cally in that the potential introduces no scales into the 
problem — its range is irrelevant since it is so short as to 
be taken to be zero while the scattering length is irrele- 
vant as it is so long as to be effectively infinite. In such 
a system the binding energy must be proportional to lu. 
We wish to find the constant of proportionality. 

Clearly, the entire effect of the short-range potential 
on the wave function for ro 3> R can be encoded as a 
boundary condition for the s-wave part of the wavefunc- 
tion at small r. One way to see this is to start with 
a finite ranged potential and consider it in the limit in 
which the strength of the potential, Vq, goes to infinity 
and the range goes to zero in such a manner as to keep the 
scattering length fixed. Since the range of the potential 
is taken to zero (to reflect the underlying short range of 
the potential compared to relevant scales) the boundary 
condition is at r = 0. Moreover the zero-range nature of 
the interaction implies that the boundary condition is in- 
sensitive to energy thus the zero energy behavior encoded 
in the scattering length is sufficient. The condition is 



u (r) r=0 
with uo(r) 



1 
a 



(4.2) 



r fdS1ip(r) 



4- 



We consider the infinite scattering length case at the 
outset, so this limiting procedure simply yields the log- 
arithmic derivative of the wave function with respect to 
r approaching zero at r = 0. Diagonalizing the Hamilto- 
nian for this system in a truncated basis and testing the 
numerical convergence of the truncation yields a binding 
energy of 



Eh — F n 



with F-* 



.60486 



(4.3) 



where the result is numerically accurate to the quoted 
precision. As advertised the binding energy is simply 
proportional to the spacing of the Landau levels of the 
free particles, w. 

Of course, an infinite scattering length comes about 
only through infinite fine tuning. In practice, even when 
the scattering length is quite large compared to the range 
of the interaction (as for nucleon-nucleon scattering) it is 
finite. This can play a significant role since the exter- 
nal magnetic field may be varied which in turn varies ro- 
Thus by changing the external field one can go from the 
region in which o > r to one in which a ~ r . It is 
important to consider the situation in which a ~ ?*o R 
since this broadens the applicability of the analysis. On 
general dimensional grounds it is apparent that if one 
uses u) to set the scale of the binding energy, then then 
binding energy can only depend on the dimensionless ra- 
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AE, 



finite — range 
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FIG. 1: The ratio of the binding energy to lo from Eq. (|4.4|l . 
F, as a function of a/ro. The solid dots are the results of 
numerical calculations. The dotted line indicates the asymp- 
totic value from Eq. (|4.3[1 while the dashed line is the weak 
field result of Eq (ETT)), 



tio a/?*o: 



E b = F(a/r )to . 



(4.4) 



We know from Eq. (|3.1[) that for small values a/ro, 
F(a/r ) -> 2a 2 /r§. Similarly we know from Eq. (|4.3j) 
that as a/ro gets large F approaches i 7 ^ ~ .60486. The 
function F can be computed numerically by diagonal- 
izing the magnetic Hamiltonian for the subject to the 
boundary condition in Eq. (|4.2p . This is shown in Fig. [TJ 

The analysis in this section is more general than the 
one in Sec. IIIII We have shown that in circumstances 
where the range of the potential is much shorter than the 
magnetic length ro one always has a bound state if the 
scattering length is attractive and if Coulomb repulsion 
can be neglected, and we have explicitly computed the 
binding energy in this situation. 

Clearly there are errors associated with approximat- 
ing the potential as zero-ranged and neglecting Coulomb 
repulsion. If one is in the unnatural regime where a ~ 
ro ^> i? one can estimate the characteristic size of these 
errors. From general power counting arguments in effec- 
tive field theory one expects that relative errors induced 
by treating the potential as having zero range are 



AE 



finite— rang 



O {R 2 /rl 



(4.5) 



The relative errors due to the neglect of Coulomb are 

a pCoulomb 

b — O(Q 2 fir a Log(R/r )) (4.6) 



decreases. The optimal value of the mag- 
netic field to minimize both errors will have the two types 
of errors as being of similar scales. Equating the two 
yields an error at the optimal field strength for the two 
approximations of AEb/Eb ~ (a/iR) 2 ^ 3 . Plugging in val- 
ues typical for the range of the potential yields an error 
of order 10% for the diproton. Thus one expects the 
analysis done here should have a regime of approximate 
validity for the diproton case. While it is possible that 
the corrections are large enough as to yield numerically 
significant corrections to the result in Fig [TJ it seems 
quite likely that a region does indeed exist for which two 
protons bind. 

The preceding error estimate was done assuming the 
scattering length is unnatural and of order ro. When the 
scattering length is natural (i.e. a ~ R) the relative size 
of the Coulomb correction is much larger and given by 
Eq. (|3.9p . Clearly the two results must smoothly inter- 
polate into one another. It is easy to understand qual- 
itatively how this comes about. If one starts with the 
extreme limit of a = oo the error estimate in Eq. (|4.6p 
remains valid since then ro sets the full length scale in 
the problem. As one decreases a with fixed ro two ef- 
fects occur: the binding energy drops as seen from Fig.Q] 
This in turn leads to a larger wave functions and hence 
smaller Coulomb corrections. However, the Coulomb en- 
ergy drops much more slowly then the binding energy. 

Ideally, one would like to calculate this in the regime 
of validity of an effective field theory treatment — i.e. 
ro 3> R, a R with ro/a arbitrary. Unfortunately, this 
can not be done without adding additional information 
about the system at short distance. If one took seriously 
our zero-range approximation for the interaction, the per- 
turbative estimate of the Coulomb contribution would 
diverge logarithmically at short distances due to the sin- 
gular nature of the Coulomb interaction. Of course, in 
nature this behavior is regulated by the short distance 
physics of the interaction. This yields a logarithmic de- 
pendence of the ratio of the short distance length associ- 
ated with the potential and the natural magnetic length 
scale. Thus, it is apparent that one needs more data 
about the interaction than merely the scattering length 
to compute the Coulomb energy. This is unfortunate 
in that one would like to test the qualitative estimates 
given in this section to ensure that there are not anoma- 
lously large numerical coefficients which might spoil the 
window where both approximations are reasonably good. 
We have done simple model calculations to explore this 
possibility. The models are consistent with our qualita- 
tive estimates; none of the models raises doubts about 
whether two protons bind in the regime where ro/a ~ .2. 



SPIN AND THE GROUND STATE 



There is a tension between the two approximations 
which is seen clearly from Eqs. (|4.5[) and (|4.6[) . While 
b e-range increases with increasing magnetic field, 



Heretofore we have neglected spin in our discussion. Of 
course, for the general case particles with non-zero spin, 
the inter-particle interactions will be spin dependent (e.g. 
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there can be a tensor force). As noted previously, for a 
generic interaction whether or not a magnetic field will 
yield binding depends on the details of the interaction 
and must be assessed on a case by case basis through 
direct simulation. As in the case of spinless particles, if 
the system is in a regime in which the Coulomb interac- 
tion can be neglected and the potential is attractive in 
all channels and all distances a bound state must occur. 
The previous variational argument holds with only mi- 
nor modification. The only new issue is that in setting 
up a variational ansatz, the state should be put into some 
particular spin state consistent with the statistics of the 
problem. 

Similarly for ro 3> R, the analysis of the previous sec- 
tions should hold for the energy of a bound state both 
for the natural scattering length case of Eq. (|3.ip and 
for the unnatural case discussed in Sect. |lVl Again, the 
only new issue is the need to consider a particular spin 
state consistent with the statistics of the problem. If for 
example one is studying identical spin 1/2 fcrmions — as 
in the case of two protons — either the spin state is sym- 
metric (spin one) and the space is antisymmetric (odd 
L) or the spin state is antisymmetric (spin zero) and the 
space symmetric (even L). Since we have been focusing 
on the case where ro R where only the s-wave inter- 
action is relevant the analysis of the previous sections 
holds provided we restrict our attention to the total spin 
zero sector. This our previous analysis shows that two 
protons with total spin zero will form a bound state in 
sufficiently strong magnetic fields. 

However, there is one critical way in which spin can 
become important in this problem. In free space it is 
always true that if the two-particle system has a bound 
state, then the ground state cannot be in the contin- 
uum. However, one cannot immediately generalize the 
intuition gleaned from this fact to the case of bound 
states in constant magnetic fields. Note that the char- 
acteristic scale of the binding energy is of order ui. The 
interaction energy of the magnetic moments of particles 
with nonzero spin with an external magnetic field are of 
the same characteristic order (assuming the g-factors are 
of order unity). Moreover, for the case of identical spin 
1/2 particles these two effects are in competition. The 
binding effect of the potential lowers the energy for the 
spin zero state; the magnetic interaction lowers the en- 
ergy for the spin 1 state (with the moments aligned with 
the external field). If this magnetic energy turns out to 
be larger then the ground state is not the bound state. 

As it happens, this is precisely the situation for the 
diproton. The interaction energy of the magnetic mo- 
ment of the two protons in a spin singlet state is, of 
course, zero. In the spin aligned state it is: 

A£ , sp i„-aiig„od = _ gBo ^ n = _9_u_ ^ hX) 

where hn the nuclear magneton is The g- factor 

for the proton is g p w 5.586 yielding A.E s P in - ali s nod w 
— 2.793 u>. In contrast, one sees from Fig. [T] that the bind- 



ing energy is bounded from above by F^lo ps .60486 uj. 
Thus, the energy the system gained by binding in the spin 
singlet channel is always more than four times smaller 
than the energy the system gains by a spin flip into the 
spin aligned total spin 1 state. 

The upshot of this is that although there is rather 
strong evidence that the diproton has a bound state in 
sufficiently strong field, the ground state will be in the 
two particle continuum. 



VI. DISCUSSION 

As noted previously, the argument that in the regime 
where Coulomb repulsion can be neglected attractive po- 
tentials for identical charge particles will lead to bound 
states in magnetic fields is quite general. However, get- 
ting into the regime in which Coulomb repulsion is neg- 
ligible requires extraordinarily large magnetic fields. We 
see from Eq. (|3.9[) that unless there is fine tuning of the 
scattering length, the Coulomb correction becomes small 

when (q 2 (/iR)^^J < 1. A quick and dirty estimate 

for systems of atomic ions (with R of order the Bohr ra- 
dius and fi ~ M p gives B$ ^> 10 11 T. Clearly this is far 
beyond the regime for which the ions are essentially un- 
altered from free space and the interaction between the 
atoms is given by its free space value. The situation only 
worsens for systems of molecular ions. 

For this reason it is apparent that if the phenomenon 
occurs in a clean system it is most likely to be for strongly 
interacting particles. The diproton system is particularly 
favorable since the scattering length is unnaturally large. 
From Eq. (|4.6[) . it is apparent that one requires magnetic 
fields of order 10 12 T for the Coulomb repulsion to be 
modest enough to be treated as a correction. This scale 
is of physical interest. Magnetic fields of this scale are 
believed to be found in nature — in particular on the sur- 
face of magnetars. Thus the novel type of bound state 
discussed in this paper appears to be physically relevant. 

One important open issue associated with this phe- 
nomenon for the diproton is whether there is a maxi- 
mum magnetic field for which the system binds. As one 
increases the magnetic field there are corrections to the 
effective field theory type treatment arising from shorter 
distance effects; assuming such corrections are modest, 
the scale of these corrections are given in Eq. (|4. 5|) . As 
Bq is increased further, however, these corrections will 
become large and the expansion will break down. In 
principle one can simply solve the problem using "re- 
alistic" nuclear potentials. However, as Bq increases, ro 
decreases. When ro becomes short enough that it is sim- 
ilar to the spatial extent of the nucleon one expects that 
magnetic effects would begin to cause substantial modi- 
fications of the nucleon substructure; this in turn could 
lead to substantial modifications of the nucleon-nucleon 
potential. Since the spatial extent of the nucleon is very 
similar to the range of the nucleon-nucleon potential it is 
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not clear whether an explicit calculation with a realistic 
potential would provide any reliable information beyond 
what can be studied via the effective theory here. 

Another important open question is whether there are 
any other strongly interacting systems besides the dipro- 
ton for which this novel class of bound state will form. As 
noted in Sect. IIIII for systems with "natural" scattering 
lengths of order R it appears that typical nuclear sys- 
tems are marginal at best. If the field is strong enough 
to overcome Coulomb repulsion, the system may well be 
out of the weak field limit. As the field increases beyond 
this limit, the magnetic interactions may alter the inter- 
nal structure of the nuclei altering the potentials. Thus, 
apart from diprotons in magnetar strength fields we know 
of no other clean examples of the phenomenon. 

In looking for other systems, it should be noted that 
this paper only considers the case for which the system 
consists of identical charged particles; this leads to a 
decoupling of the center of mass and relative motion. 
While this simplifies the problem both conceptually and 
computationally, it also forces us to consider very strong 
magnetic fields to overcome the effects of the Coulomb 
repulsion. One might wish to consider the situation in 
which a charged particle and a neutral particle interact 
in the presence of a magnetic field. In such a situation, 
Coulomb repulsion would no longer be an issue. How- 
ever, the notion of a two particle bound state becomes 
a bit more involved in the presence of a magnetic field 
since the relative and center-of-mass motion are coupled. 
In practice, a useful definition is that the two particles 
are in an energy eigenstate with the property that the 



position of the two particles are strongly correlated — if 
one observes one of the particles at some point, then the 
probability of not finding the other particle a distance r 
away decreases exponentially with r for large r in any 
direction. An alternative definition is that there exists 
a state with an energy less than the lowest Landau level 
for the isolated charged particle. The simple qualitative 
arguments as to why one expects such bound states a 
priori are weakened as compared to the case of identi- 
cal charged particles; we do not know whether magnetic 
fields can cause a charged and neutral particle acting 
through a potential to bind (in the sense defined above) 
when they do not bind in free space. We will consider 
this possibility in a subsequent work. 

Finally, we turn to the implications of this result from 
the astrophysical perspective. On the one hand, the 
magnetic field strength in magnetars appears to be large 
enough for bound diprotons to form. An optimistic view 
on this would be that this can have important effects 
on the dynamics of the magnetar leading to some obser- 
vational signature. On the other hand, as discussed in 
Sect. [V] the bound diprotons arc not the ground state 
of the system; spin aligned unbound protons are. To 
the extent that protons will tend to stay in their ground 
state, this suggests that it will be difficult to find an ob- 
servational signature of bound diprotons. Our results, 
however, suggest that nuclei that are unstable due to 
excessive number of protons may become stable in the 
presence of a magnetic field of this magnitude. This shift 
of the nuclear valley of stability towards higher values of 
Z deserves further study. 
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